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Abstract 

We compute the two-loop term in the perturbation series for the quark-mass in 
the lattice Heavy Quark Effective Theory (HQET). This is an ingredient in the 
matching factor required to obtain the 6-quark mass from lattice simulations 
of the HQET. Combining our calculation with numerical results from the APE 
collaboration, we find, at two-loop order, m;, = m^ LS (m^ ls ) = 4.41 ± 0.05 ± 
0.10 GeV. It was expected that the two-loop term would have a significant 
effect, and this is indeed what we find. Depending on the choice of "reasonable" 
coupling constant in the one-loop estimates, the result for can change by 
several hundred MeV when the two-loop terms are included. 
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1 Introduction 



In this paper we evaluate the two-loop perturbative matching factor required to obtain 
the mass of a heavy quark from lattice simulations of the Heavy Quark Effective Theory 
(i.e. from simulations of static heavy quarks). This is the next-to- next-to-leading-order 
(NNLO) term in the matching factor. We combine this matching factor with numerical 
results from simulations to obtain: 



m b = 4.41 ± 0.05 ± 0.10 GeV 



where 



(1) 



(2) 



at NNLO. The first error in eq. (0) is due to the uncertainties in the values of the lattice 
spacing, in the value of the strong coupling constant and in the numerical evaluation of 
the functional integral. The second error is an estimate of the uncertainties due to our 
ignorance of 3-loop and higher order perturbative terms. 
For the Lagrangian of the HQET we take: 



•-HQET 



hD. 



i + r 



h 



(3) 



where h is the field of the heavy quark, and we use the following definition of the covariant 
derivative D: 



h(x) = U^(x)h(x +{/,) — h{x) 



(4) 



where x + ft is the neighbouring point to x in the /i-direction. U^(x) is the link variable 
from x to x + fi. For the light-quark action we take an "improved" generalisation of the 
Wilson action. For each quark flavour the action is: 



S q — Sw + csw 5sw 



(5) 



where Sw is the Wilson fermion action: 



S * = Y, \ ~ o E k 3 ^ 1 - jMxWx +A) + 1>(x + A)(l + l,)Ul(x^(x) 



+ ip(x)(m + 4)ip(x) 



(6) 



and Ssw is the Sheihkoleslami-Wohlert (or "clover") action || 

Ssw = ~ ^2 \^{x)a^ u F^(x)il)(x) 



(7) 
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(a) (b) 

Figure 1: (a) A one- loop diagram and (b) a two- loop diagram contributing to 5m. The 
solid lines represent static heavy quarks. 

F^y is the lattice expression of the field-strength tensor obtained by averaging the four 
plaquettes lying in the (/i, v) plane and stemming from the point x. The coefficient csw is 
equal to 1 at tree-level (with csw — 1 an d using appropriate operators the discretisation 
errors in physical quantities are reduced from 0(a) for the Wilson action to 0(a s a)). 
Recently much effort has been devoted to fixing csw non-perturbatively in such a way that 
the errors are reduced to ones of 0(a 2 ) 0. For the perturbative calculations in this paper 
we will keep csw and Nf, the number of light-quark flavours as free parameters Q. 

The starting point for our evaluation of the heavy quark mass is the computation of 
the correlation functions of the time-component of two axial currents: 

C{t) = X>|A 4 (x) 4(0)|0> (8) 

~ Z 2 exp(-£t) , (9) 

where we have assumed that t is positive and sufficiently large so that we can neglect the 
contributions from excited states. The current takes the form A±(x) = h(x)^4^^q(x) (the 
factor of 74 can be replaced by the Identity since it is adjacent to the heavy quark field). 
By fitting the correlator C(t) as a function of t, both the prefactor Z and the exponent 
8 can be determined. From Z we can obtain the decay constant of a heavy pseudoscalar 
meson in the static limit, whereas from £ we obtain m& ||. Matching QCD and the Lattice 
HQET one finds that §: 

m P° le = M B - £ + 5m , (10) 

where m^ ole is the "pole" mass of the tVquark, Mb is the physical mass of the S-meson 
and 5m is the mass generated in the static theory, eq. @, in perturbation theory (a5m = 
a s X + a 2 Xi, where the coefficient X\ depends on the choice of coupling constant, a s , 
used in the perturbative expansion). Even though there is no explicit mass-term in the 
bare action of eq. (§), Feynman diagrams, such as those in fig. [I], generate such a term 
which is proportional to 1/a, where a is the lattice spacing. Once ml° le is known, then 

1 The results are also presented for an arbitrary number of coulours, N. 

2 In practice it is generally advantageous to use "smeared" interpolating operators in order to enhance 
the overlap with the ground state. 
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any physical, short- distance, quark mass can be determined using continuum perturbation 
theory, e.g. 



m b = (M B -£ + 5m) 



1 _ | - („. 66 - 1M N f ) (^-J + O (a?(m t) ) 

(11) 

where Nf is the number of light-quark flavours and a s is defined in the MS renormalization 



scheme. The term in square brackets in eq. (|TT|) is the pertubative expansion of the factor 



which relates the pole mass to the MS one ||] . The main aim of this paper is to calculate 
5m to two-loop order. Using the MS coupling a s (m b ) as the expansion parameter we find 

a5m = 2.1173 a s (m b ) + { (3.707 - 0.225 iV/) \n(m b a) 

- 1.306 — Nf (0.104 + 0.100 c sw — 0.403 Cg W ) } ot 2 s (rni) ■ (12) 

Consider eqs. (|10"D, ( [TTj) and (|12|). Neither the perturbation series for 5m nor that 
relating m^ ole and m b in the square brackets in eq. (|TT|) is Borel summable and both 
contain renormalon ambiguities. In each case the leading ambiguity is of O(Aqcd)- Since 
m£ ole is not a physical quantity, the presence of a renormalon on the right-hand side of 
eq. (|T0"D is not inconsistent. m b , however, is a physical quantity, and the leading renormalon 



ambiguity cancels when the two perturbation series in eq. ( [TT]) are combined ||, || |7j f\. 
The remaining renormalon ambiguity is of 0(AQ CD /m b ) which is beyond the precision we 
are considering in this paper. Related to the problem of renormalons is that of power 
divergences (i.e. the terms which behave as powers of a -1 ). 5m diverges linearly with the 
inverse lattice spacing, cancelling the linear divergence in 8 (this cancellation is necessarily 
only partial since we truncate the perturbation series for 5m at a finite order of perturbation 
theory). Thus by computing £ numerically, and calculating 5m and the relation between 
the pole-mass and m b in perturbation theory (this relation is known to two- loop order H), 
m b can be obtained. It is free of power divergences and renormalon ambiguities. 

Although we do not present the discussion explicitly in the formalism of an Operator 
Product Expansion, the calculation of the quark mass here is an example of the evaluation 
of a power correction. To lowest order in the heavy quark expansion we can take m b = Mb, 
and we are evaluating the first power corrections to this (i.e. the correction of O(Aqcd) 
relative to the leading term). In ref. we argued that the evaluation of power corrections 
is difficult, requiring high-orders of perturbation theory in order to subtract the power di- 
vergences to sufficient precision. To illustrate this point for the heavy-quark mass, we show 
in subsection |3.1.3| that the one-loop result changes by several hundred MeV depending on 



which "reasonable" coupling constant is used in the one-loop result (specifically we use the 
MS coupling at an average momentum q*, the V coupling defined from the potential also at 
the scale q* and a boosted lattice coupling constant). The two- loop calculation described 
in this paper reduces this uncertainty to 100 MeV or less. 

3 In order for this cancellation to be manifest, a single coupling constant for both perturbation series 
must be used (we use a s (rnb)). 
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In the present calculation the ultra-violet divergences are linear in the inverse lattice 
spacing. For other physical quantities for which the power divergences are quadratic (or 
even higher order) it is even more difficult to control the perturbative corrections |7|. 

The plan for the remainder of this paper is as follows. In the next section we present 
the outline of the perturbative calculation of 8m up to two-loop order. In an attempt 
to make the steps of the calculation clearer, we have relegated the technical details to a 
(rather long) appendix. In section ||] we use our perturbative results in order to determine 
the mass of the 6-quark from values of S obtained in numerical simulations. Finally in 
section || we present our conclusions. 



2 Calculation of 5m to Two-Loop Order 

In this section we outline the evaluation of 8m up to two-loop order in perturbation theory. 
Instead of evaluating self-energy diagrams, such as those in fig. [1], directly, we exploit the 
fact that two-loop perturbative results for rectangular Wilson loops have already been 
presented ||, ||. We therefore determine 8m from the exponential behaviour of large 
Wilson loops as we now explain. 

Consider a Wilson loop (W(R,T)) in the \x-v plane, of length R and T in the \x and v 
directions respectively. Since, in perturbation theory, the potential between static-quarks 
falls like 1/r with their separation r, the expectation value of large Wilson loops decrease 
exponentially with the perimeter of the loops, specifically 

( W(R,T) ) ~ exp(-28m(R + T)) . (13) 

Hence the perimeter term in log(( W(R, T) )) is simply —28m(R + T). 

Following ref. || we define the first two coefficients of the perturbative expansion for 
(W(R,T)) as follows: 

( W(R, T) ) = 1 - g 2 W 2 (R, T) - g 4 W 4 (R, T) + 0(g 6 ) , (14) 



where g is the bare lattice coupling constant. By using eq. fll3|) we obtain the perturbative 
expansion for 8m in terms of these coefficients, 



8m 



2(R + T) 



g 2 W 2 (R,T) + g\W 4 (R,T) + 1 1 Wi(R,T)) . (15) 



We now evaluate W 2 and W4, only keeping terms which grow at least linearly with the 
dimensions of the Wilson loop. 

2.1 Evaluation of W 2 (R, T) 

In ref. || we find the following integral representation of W 2 for an S'L r (iV)-gauge theory: 
W (R T) - {N2 ~ 1} T ^ ^(iP^sm^R) f 1 1 \ 
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where 



D(p)=4j>n 2 (ip A ). (17) 



A=l 



Consider the first term in the integrand in eq. (16), i.e. the term containing a factor of 
1/ sin 2 (|p Ai ). By inspection (i.e. by using power counting) it can be seen that the region 
of small p^ gives a contribution proportional to R (and similarly the second term gives a 
contribution proportional to T). By changing variables from p^ to z = exp(ip^) (so that 
the contour of integration in z is the unit circle), and studying the singularities of the 
integrand in z, the integration over z can be performed. In this way one obtains 

(N 2 -l)— (N 2 -l) r 
W 2 (R, T) ee { -^^W 2 (R, T) = { -^ r 1 {(R + T) X - Y{R) - Y(T)} , (18) 



where X and Y are three-dimensional integrals defined in eqs. (|A.1|) and ( |A.4|) of the 



appendix. X is a finite integral, whereas Y(R) grows logarithmically with R. 5m at 
one-loop level is given by the term in the perimeter which is proportional to X. We will, 
however, require the terms proportional to Y(R)+Y(T) when we evaluate 5m at two-loops 
(see eq. ((L5|) ). Equation fll8|) implies that 

AT 2 - I 

Sm^-^Xg 2 + O(g'), (19) 
which is the well known one-loop result. 

2.2 Evaluation of W 4 (R, T) 

The evaluation of Wi{R,T) is the main calculation of this paper, and we describe this 
calculation in some detail in the Appendix. We write 

W 4 {R, T) = Wf non (R, T) + Wt mion (R, T) , (20) 

where Wf non (R,T) is the contribution from diagrams in the pure-gauge theory and 
Wl eimion (R,T) is that from dia grams containing light-quark loops. The purely gluonic 
contribution is given in eq. ( |A.50|) 



Wt°\R, T) = - {N l 2N l ? ' X\R + Tf + {N w l ? X (Y(R) + Y(T)) (R + T) 

(N 2 - 1) 

+ 1 ; {3X 2 - 2XL + 12W} (R + T) 

(21) 
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1 TITOO'TQ 1 

llltcg,! cLl 


.L/dllltiU. Ill t:L^. 


V cLlllt: 


x 


el a ih 




L 


(M) 


0.30987 


W 


(M) 


1.313(3) • 10~ 2 


V x 




7.008(5) • 10~ 2 


v 2 




-2.130(1) • 10~ 2 


V w 


dA.14j) 


-2.145(1) • 10~ 2 


T/a 


OA.isp 


2.63(2) ■ 10- 3 


V b 


(IA.15I) 


-2.974(1) ■ 10~ 2 



Table 1: Values of the integrals used in the evaluation of 5m to two- loops in perturbation 
theory. 



and the contribution from graphs containing light-quark loops is given in eq. ( A.52Q 



W { ~(R, T) = {N \ 6 ^ )Nf (R + T) (Vw + V m ) , (22) 

where Nf is the number of active light-quark flavours. The integrals X, L, Y, W, Vip, Vw 
and Vsw m these expressions are defined and evaluated in section |A.1| of the appendix. 



Vsw is the contribution from the improvement term in the fermionic action and we write it 
as Vsw = Vsw c sw + V s b w c| w . With the exception of Y(R), which depends logarithmically 
on R and cancels out of the expression for 5m, we present the numerical values of these 
integrals in table [l]. 



2.3 The Perturbation Series for 5m 

Having determined the relevant terms in Wv{R, T) and W^(R, T) we can readily obtain 
the perturbation series for 5m using eq. (pl|). At one- loop order, the result is presented in 
eq. ( |i~9"D which, setting the number of colours N equal to 3, we rewrite as: 

a5m ~ 2.1173a + O(a 2 ) , (23) 

where the subscript on a reminds us that the expansion is in terms of the bare coupling 
constant in the lattice theory. 

The two-loop contribution to 5m is proportional to W^(R, T) + |PV 2 2 (i?, T) (see eq. ( |T5D ): 

Wi(R, T) + \W*{R, T) = (R + T) ^ ~ ^ I 3X 2 + 12W + 6Vi + 12 V 2 4 

192 

3(iV 2 -4) l2N ft 
12LX+ y m2 ) X + ^(V W + Vsw)} • (24) 
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The terms quadratic in R and T explicitly cancel in the combination W4 + \W%, as do 
those proportional to the integrals Y(R) and Y(T) as they must. This is an important 
consistency check of our calculations. All the integrals on the right-hand side of eq. ( p4"|) 
are finite and independent of R and T. Substituting N = 3 and the numerical values of 
the integrals (see the appendix) in the right-hand side of eq. (0), we obtain f\ 

Wi(R, T) + §Wf T) = (R + T) {0.14124 + ^(-0.00358 + 0.00044 c sw - 0.00496 c 2 sw )}. 

(25) 

We thus arrive at the two-loop expression for 5m: 



, 5m ~ 2.1173 a + {11.152 + ^(-0.282 + 0.035 c sw - 0.391 c§ w )} ag + O(ajj) . (26) 



Eq. (|26| ) is the main result of this paper. In the following section we will exploit this result 
to extract the mass of the 6-quark from the values of £ computed in lattice simulations. 



3 Lattice Computation of nib 



In this section we discuss the evaluation of the heavy quark mass, from lattice computations 
of £ using the two- loop result for 5m in eq. (|2HD . For illustration we will determine the 
MS-mass, m&, defined in eq. (0), from which any other short-distance definition of the 
heavy-quark mass can be obtained using continuum perturbation theory. 
As discussed in the introduction, the relation between m& and £ is: 



m b = (Mb — £ + 5m) 



_ 4 / a s {m b ) 
3 V vr 



- (11.66 - 1.04N 



f) 



7T 



(27) 



where Nf is the number of light-quark flavours and a s is defined in the MS renormalization 
scheme. In eq. (0), Mb = 5.279 GeV is the physical mass of the 5-meson, £ is computed 
in lattice simulations and the remaining terms are calculated in perturbation theory. The 
term in square brackets in eq. (|ll|) is the perturbation expansion of the factor which relates 
the pole mass to the MS one 0. 

The two-loop expression for 5m in eq. ( p6|) is given in terms of «o, the bare coupling 
constant in the lattice theory. In order to achieve the explicit cancellation of renormalon 
singularities the same coupling constant needs to be used in the expansion of 5m and in 
the relation between the pole mass and the MS mass (i.e. in the expression in square 
parentheses in eq. (|li"D ). The relation between the MS coupling a s and the lattice coupling 
a is: 



n.„ I - J = a + di(s) a + d 2 (s) a + 



(2f 



4 The errors in the coefficients in eq. (^5|) (and following equations) are due to those in the numerical 
values of the integrals given in the appendix. They are typically 1 (or possibly 2) in the last digit, and are 
negligible in the evaluation of . 
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where 



dAs) 



^- + 2.13573iV + N f (-0.08413 + 0.0634 c sw - 0.3750 c 2 sw ) 



and fa = (UN - 2N f )/3. 
gauge theory without fermions is found in ref. 
fermions (i.e. with csw — 0) can be obtained from ref. JTT 



(29) 

The contribution to the right-hand side of eq. (|29"D for the 
10] , and the quark contribution for Wilson 
The coefficients of csw and c| w 



are presented here for the first time. Although not needed for our present calculation, the 
two- loop coefficient in eq. (|28|), d 2 (s), for the pure-gauge theory can be found in ref. |12| . 

Combining eqs. (2^) and (|2~9"D we obtain the result in eq. fll2T) for 8m expressed in terms 
of the MS coupling. 



3.1 Numerical Results 

The results for £ computed in lattice simulations to date, have been obtained in the 
quenched approximation. For the purposes of this study we take the results from the APE 
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a£ = 0.61(1) at /3 = 6.0 (a" 1 = 2.0 (2) GeV) (30) 
a£ = 0.52(1) at = 6.2 (a" 1 = 2.9 (3) GeV) (31) 
a£ = 0.460(7) at /3 = 6.4 (a" 1 = 3.8 (3) GeV) . (32) 

Of course the quenched calculation is incomplete and there is no procedure for determining 
rrib from quenched computations of £ which is totally satisfactory. We will now describe 
our approach, but note that many readers may choose to follow different procedures, which 
are equally valid in the absence of unquenched results. We take the quenched results in 
eqs. (|30| ) -(|32| ) and determine the pole mass by combining them with the perturbative 
result for 8m written in terms of the MS coupling constant at the scale % obtained using 
eqs. ( pgj) and (p9"D with Nf = 0. We then derive m^ from this value of the pole mass using 
continuum perturbation theory, eq. (|TTD, but also with Nf = 0, to ensure the cancellation 
of renormalon singularities. We present the details and results in subsection |3.1.1| . 

In subsection [3.1.2| below we discuss the fermionic contributions in perturbation theory 



and demonstrate that they are large (of (9(100) MeV or more for the lattice spacing con- 
sidered in this paper). Finally in this section we compare our two-loop results with some 
standard procedures used to try to optimise the one-loop results, such as using the boosted 
coupling constant or a coupling defined at some average momentum for the process, q*, 
(see subsection |3.1.3j ). 



3.1.1 Quenched Results for m b 

In the quenched approximation we can rewrite 5m in the form 

a5m = 2.U73a s {fn b ) + (3.7071n(m fe a) - 1.306) a 2 s (m b ) . (33) 
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p 


a' 1 (GeV) 


m b (GeV) 
eq. (0) 


m b (GeV) 
eq. © 


6.0 


2.0(2) 


4.37(5) 


4.45(5) 


6.2 


2.9(3) 


4.38(7) 


4.47(7) 


6.4 


3.8(3) 


4.38(7) 


4.47(8) 



Table 2: Values of m b obtained using eqs. (0) and fl35l) . 



Most of the very large coefficient 11.152 in the perturbation series written in terms of «o 
(see eq. fl2lf) ) has been absorbed into the coupling ot s . This is frequently the case JjlJ]. The 
implicit equation for m& is (see eq. ([TT|) ) 



rn.{, = A 



1 + — {2.1173 a s (m b ) + (3.707 ln(m 6 o) - 1.306) c^(m b )j 
Aa J 



x 



1 - 



Aa s (m b ) 



11.66 



7T 



7T 



(34) 



where A = tub — £■ To the order we are working in, eq. (|34|) is equivalent to: 



m b = A 



2 1173 

1 + < — : 0.4244 } a s (m b ) 



Aa 



+ { -T- (3.707 ln(m 6 a) - 2.207) - 1.181 j- a 2 Jm b ) + ... 
Aa 



■ (35) 



and 



will be an indication of the 



The difference of the results obtained using eqs. 
error due to 3-loop and higher order terms. 

We estimate m b using equations (0) and (^). At this stage, in spite of using the 
quenched results for S, we simply assume that the result is the physical one and for 
a s (jn b ) we take the coupling constant obtained from a s (M z ) = 0.118(3) using the two- 
loop evolution equation with five flavours (this gives a s (m b ) ~ 0.22(1)). Eqs. (|34"D and 
(j3~5|), which contain m b both on the left and right-hand sides, are solved by iteration and 
the corresponding results are presented in table 0. As mentioned above, the difference of 
the results obtained using eqs. ( p4|) and (|35| ) are due to our ignorance of three-loop and 
higher-order perturbative corrections. This difference is about 80-90 MeV. The results for 
different values of the lattice spacing are in remarkably good agreement. The errors on the 
values of fn b in table |2] are due to the uncertainties in the values of the lattice spacing, in 
the values of S and in a s (Mz)- The dominant component in the errors is the first one, and 
is a consequence of the fact that using different physical quantities to determine the lattice 
spacing in quenched simulations, leads to differences of O(10%) or so in the estimates of a 
(which correspond to differences of about 50 MeV in fn b ). The errors in the results for m b 
in table are strongly correlated and should not be combined in quadrature. For our best 
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p 


a' 1 (GeV) 


m b (GeV) 


m b (GeV) 


eq. (@) 


eq. © 


6.0 


2.0(2) 


4.52(4) 


4.60(4) 


6.2 


2.9(3) 


4.63(5) 


4.75(5) 


6.4 


3.8(3) 


4.78(6) 


4.92(7) 



Table 3: Values of m b obtained using eqs. (|34"D and (|35"D in each case dropping the terms 
which are explicitly of 0(a 2 s (m b ) ). 



value we take the results at (3 = 6.0 and quote: 

m b = 4.41 ±0.05 ±0.10 GeV . (36) 

The first error in eq. (^) is just that in table and the second is an estimate of the 
uncertainty in higher order perturbative corrections, as manifested in the different values 
in the third and fourth columns of the table. 



3.1.2 The Unquenched Calculation 

We now briefly comment on the result including fermion loops. Consider eq. (JT2|). In 
this case it is not true that most of the fermionic two-loop contribution is reabsorbed into 
the MS coupling constant. The two-loop fermionic contribution is large, particularly with 
the improved action. Even with the tree-level improved action (csw — 1) the fermionic 
two-loop contribution to 5m varies from about 80 to about 150 MeV for lattice spacings 
in the range considered in this paper. For full 0(a) improvement one expects csw> and 
hence the two-loop contributions, to be larger still. Since the fermionic contribution to 5m 
is large, it is likely that unquenched results for S will be significantly different from those 
in eqs. (j30l)- (|32|) for the corresponding values of the lattice spacing. Of course we are not 
yet in a position to check whether this expectation is true. 



3.1.3 Comparison with the One-Loop Result for nib 

In table |] we present the results for m b obtained using one-loop perturbation theory in 
terms of a s (m b ). The errors in table |3] are also those due to the uncertainties in the values 
of the lattice spacings, £ and a s (M z ). Using this procedure, the variation of the results 
with (3 is relatively large and it is difficult to extract a meaningful result for m b . 

In many calculations in which the perturbative matching factors are known only to 
one-loop order a different procedure is followed. Instead of using a s (m b ) as the coupling 
constant in lattice perturbation theory, one uses a s (q*) where q* is an estimate of the 
"typical" loop momentum in the process |I3 |. We choose to define q* by: 



'°si< a «*> i = x J ^mr • (37) 



n 



where X and A(k) are defined in eqs. ( |A.1| ) and ( A.2j ) respectively. With this definition 



q*a = 1.446. Alternatively one might use a "boosted" lattice coupling constant [14]. In 



order to get an estimate of the precision and reliability of such procedures we define 

e(5m) = (5m) two _i oops - (5m) one _ioo P , (38) 

where (<5m) two -ioo P s is the two-loop result given in eq. ( |33| ) and (5m) one _i 00 p is the one-loop 
result obtained using one of the procedures described above. We find the following: 



i) The one-loop result 5m = 2.1173 a s (q*) a -1 reproduces the two-loop result in eq. (031) 
remarkably well when the MS coupling constant is used. Specifically e(5m) ~ 10, 15 
and 35MeV for (3 =6.0, 6.2 and 6.4 respectively. 

ii) The one-loop result varies significantly with the choice of "reasonable" coupling. For 
example, using the coupling defined from the inter-quark potential, ay(g*), one finds 
e(5m) ~ -210, -235 and -315 MeV for (3 =6.0, 6.2 and 6.4 respectively Q. 

iii) The use of a boosted coupling constant, a s = ao/uQ, where ao is the lattice coupling 
and uq is some estimate of the average value of the link- variable, underestimates the 
perturbative corrections significantly. For example, using the trace of the plaquette 
to define u , e(5m) ~ 510, 610 and 670 MeV for (3 =6.0, 6.2 and 6.4 respectively. 



This is the procedure which was used in ref. |T5|], and is the main reason for the 
low central quoted in that paper (mb = 4.15 ± 0.05 ± 0.20 GeV), where the error of 
200 MeV was an estimate of the effects of higher orders in perturbation theory. 

Thus we see that the range of one-loop perturbative results for 5m, and hence for m~b is 
very large; indeed it is considerably larger than the precision required for the lattice results 
for the quark mass to be phenomenologically interesting. This is the main motivation for 
the two-loop calculation presented in this paper. 

It should also be noted that the two-loop result in the (continuum) perturbation series 
relating mf ole and m^ is also large (changing the result for by 0(250 MeV) ). Resum- 
mation techniques can also be used to try to improve the convergence of this series (see 
for example refs. |16[ and references therein). 



4 Conclusions 

The mass of the 6-quark is one of the fundamental parameters of the standard model of 
particle physics. It can be determined from lattice simulations of the HQET, by computing 
the time behaviour of hadronic correlation functions. An ingredient of such a determination 
of Mb, is the perturbation series for 5m, the mass generated in the lattice formulation of 
the HQET. In this paper we have evaluated the two-loop term in the lattice perturbation 

5 For the numerical value of ay(q*) we take ay{q*) — a s (q*)(l + 0.822 a s (q*) ), and use the physical 
value of the MS coupling a s (q*) obtained by evolving the result a s (Mz) = 0.118(3) using the two-loop 
renormalisation group equation. 
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series for 8m, and the result is presented in eq. Combining our result for 8m with 

the values of £ obtained by the APE collaboration, we obtain the result for given in 
eq. ©. 

We would like to determine with a precision better than O(Aqcd)- This is an 
example of a calculation of a power correction (at leading order we can use the mass of the 
pseudoscalar meson to estimate the quark mass). As with any power correction, to achieve 
the required precision is very difficult 0]. We have argued in section |3.1.3| that, using 



one-loop matching, the uncertainty in the determination of nib is several hundred MeV. 
Our estimate (which in reality can only be an educated "guesstimate" ) of the uncertainty 
in nib after the two- loop effects are included is that it is of 0(100 MeV). 

We consider the question of the precision which can be reached when only using one- 



loop matching to be so important that we repeat here part of the discussion of section p. 1.3 
At one-loop level we have 

8m = 2.1173 a s a' 1 . (39) 

As an example let us consider the simulation at (3 = 6.0 (for which a -1 = 2.0(2) GeV). 
Now as "reasonable" choices of a s we take 

i) the MS coupling at q*, a s (q*) ~ 0.253, so that 8m ~ 1.07 GeV. 

ii) the coupling in the "potential" -scheme at q*, ay (g*) = 0.306, so that 8m ~ 1.29 GeV. 

iii) the "boosted" coupling, defined by a /iiQ ~ 0.134, where u is an estimate of the 
avearge link- variable, obtained from the fourth root of the plaquette. With this 
boosted coupling 8m ~ 0.57 GeV. 

Since the perturbative coefficient is relatively large (2.1173 rather than 1/ir say) and a -1 ^> 
Aqcd, the spread of results obtained using reasonable choices for the expansion parameter 
a s is several hundred MeV. It is, therefore, not possible to achieve a precision in 8m (and 
hence in nib) which is better than O(Aqcd) without calculating the two-loop (or even higher 
order) matching coefficients. This is an example of a generic problem in the evaluation of 
power corrections, and is also not restricted to lattice computations. 

It should be said that our view that the ignorance of higher order perturbative co- 
efficients implies that the uncertainties in the results for power corrections to physical 
quantities are large is not universally accepted. In this paper we have confirmed our view 
with a specific two-loop calculation. For mb the leading correction is linear, and corre- 
spondingly we have had to subtract the linear divergence (i.e. the terms which diverged 
linearly in 1/a) in £ using perturbation theory. In other important examples one needs 
to subtract quadratic or even higher order divergences and the difficulty to achieve the 
required precision increases enormously. 

6 In the following examples we take the central value a -1 = 2 GeV. There is a 10% error associated with 
the uncertainty in the scale, but as we wish to study the variation with the choice of coupling constant, 
we keep a fixed value of a. 
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A Appendix: Two-Loop Contribution to Large Wil- 
son Loops 

Our calculation of 5m is based on the fact that it is proportional to the logarithm of the 
perimeter term in the perturbative expansion of large Wilson loops. In this appendix we 
present the details of the evaluation of the two-loop contribution to large Wilson loops. 

The appendix is organised as follows. In the opening section we define the multidi- 
mensional integrals which appear in our result for W4 and present their numerical values. 
A description of the evaluation of the terms in W4 which grow at least linearly with R 
and T in the pure gauge-theory (i.e. in the theory without light quarks) is presented in 
section A.2 . The fermionic contributions to W4 are evaluated in section A.3| . 



A.l Integrals 

The two-loop contribution to 5m will be presented in terms of a number of multidimensional 
integrals, X, L, Y, W, Vi )2 , Vw an d Vsw, which have to be evaluated numerically. In this 
subsection we define these integrals and present their values. All the integrals have limits 
of integration from -n to 7r for each momentum component. 

1. The first three integrals are one- loop ones. In each case it is straightforward to 
perform the integral over one component of momentum analytically, leaving a three- 
dimensional integral to be performed numerically. The first of these is 

/d 3 k 1 
Ml) = °- 50546 ^ 

where 

3 

A(k) = ^(1 - cos^)) . (A.2) 
1=1 

X is the integral which contibutes to the one-loop component of 5m. 

2. The second one-loop integral is 

/d 3 k 1 
7- ,„ . = 0.30987 . (A.3) 

(2tt)3 ^/{l + A{k)f-l) 

3. The third integral arises in the calculation of the wave-function renormalisation at 
one-loop level: 

y(r)B /(0 V(iVi (1 - n (A - 4) 



15 



where 



(3 = 1 + A- ^(l + Af-l , (A.5) 

and A and j3 are implicitly functions of k. Without the f3 T term on the right-hand 
side of eq. ([A.4Q the integral would be divergent (it depends logarithmically on T). 



In the final result for 5m, any dependence on Y(T) or Y(R) must cancel. 

4. The remaining integrals are 7-dimensional and for these the results are generally 
obtained with poorer relative precision. The first of these is 



W = 



d 3 k 1 f d 3 p 



(2n) 3 A(k) J (2tt) 3 A(p) 



dq sin (q) 



2?r (4sin 2 (g/2) + 2A(k))(4 sin 2 (g/2) + 2A(p)) 



1.313(3) • 10~ * . (A.6) 



5. There are two integrals which arise from the gluonic contribution to the vacuum- 
polarisation. The first of these is 

d 3 p d 4 k (l + cos(Ag) f4-C A (2p + k) 



Vl ~ J (27r)*(2n)* A*(p)(4-C 4 (k)){ 4- C±(p + k) 1 

= 7.008(5) • 10" 2 , (A.7) 

where p^ = (/i is a fixed direction and the integral over p is over the 3 components 
other than the /z-component). The function C4 is defined by 

4 

C 4 (r) = ^cos(r A ) . (A.8) 

A=l 

6. The second integral arising from the gluonic vacuum-polarisation graphs is 
d 3 p d A k sin 2 (£g / 1 + Ei^cosfe) 2 



(27r) 3 (27r) 4 A 2 (p)(4-C 4 (A;)) t 2 (4 — C 4 (p + k) ) 4 - C 4 (k) 
2.130(1) • 10~ 2 , (A.9) 



where, again, fi is a fixed direction and = 0. 

7. The final two integrals correspond to graphs containing fermionic contributions to 
the vacuum-polarisation. For Wilson-fermions the relevant integral is 



d 3 p d A k 1 1 
(27r) 3 (27r) 4 A 2 (p) S(k) 



Zw(p, k) _ Z w (p = Q,fc) 
S(q) S(k) 



(A.10) 
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where q = p + k, p^ = 0, 

4 

S(k) = Y,*™ 2 (k x ) + W 2 (k) , (A.ll) 



A=l 



W(fc)=2^sin 2 ^ , (A.12) 

A=l ^ ' 



and 



Z w (p,k) = — 8 cos 2 (/c M ) sin 2 (/c M ) + 4 sin(/c A ) sin( 



A=l 



+4 cos{2k fl )W{k)W{q) - 8 cos(fc M ) sin 2 (^) [W(fc) + W(g)] . (A.13) 



the Wilson term in the fermion propagator defined in eq. flOp , should not 



be confused with the integral W defined in eq. ( |A.6| ). We find that 

V w = -2.145(1) • 1(T 2 . (A.14) 

8. Finally we present the integral, Vsw required to evaluate the additional fermionic 
contributions when the SW-improved action is used, 

f d 3 p d 4 k 1 Z sw 
Vsw -J (2tt)3 (2tt)M 2 (p) S(kjW) ' { } 

where, again, q = p + k, p^ = and 

Zsw = c S w I 4sin 2 (fc M ) ^sin(p A ) [sin(g A ) - sin(/c A )] 

I A=l 

/ 4 4 

-4 cos(A; M ) I W(k) ^ sin(p A ) sin(g A ) - W(q) ^ sin(p A ) sin(fc A 

\ A=l A=l 

-4 W I ^2sin 2 (Ag - ^ sin(g A ) sin(A; A ) + W(k)W(q)j $>n 2 (p p ) 

4 4 "| 

+2 sin(p A ) sin(feA) ^ sin(p CT ) sin(g g ) > . (A. 16) 

A = l (7=1 J 

Evaluating the integrals numerically we find 

VW = 2.63(2) • 1(T 3 c sw - 2.974(1) • 1(T 2 c 2 w . (A.17) 
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A.2 The Gluonic Contribution to W 4 (R,T) 

In this section we outline the evaluation of the gluonic contribution to the two-loop term in 
the Wilson loop, — g A W^R, T). Following ref. |§, we distinguish 5 contributions to W4 Q: 

W 4 (R, T) = W Sl + Wj + W n + W VP + W V p . (A.18) 

We now evaluate each of these in turn, picking up the contributions which grow as R, T, 
R 2 , T 2 or RT. 



Wsv This contribution comes from the "spider" graph and is given in eq. (3.6) of ref. 
It is proportional to 



1 



d A k d A p 

(2tt) 4 (2vr) 4 D(p)D(k)D(p + k) 

sm{\p v T) sin(|(2/c + p) 



— —j — — smUp v T) cos( 5 p„) sm( 5 (2A; + p)^) 



x 



sm{U2k + p) u T) 



sm(~p„) sin(|/c v ) sin(|(p + fc)„) sin(ifc„) sin(|(p + 
sin(|p M -R) sm^k^R) sin(|(p + fc) y T) , 
smd^J smd^J sm(|(p+ 

x cos(|(p + kj^R) (xn(\(p + k)Jam(\((k-p) v ))\ + {(jJi,R) <-> (z/,T)} 



(A.19) 



where -D(p) is defined in eq. (|17|). By detailed inspection of the infrared behaviour of the 
integral, one can see that there is insufficient enhancement from the denominators to give 
a term proportional to R or T. Thus there is no contribution to 5m from Ws 1 - 



Wi: W 7 / is defined to be the "expectation value of the non-abelian part of the order 
g A term in the expansion of the Wilson loop" ||. It is presented in eq. (3.7) of ref. |§ 
and contains several terms which we consider in turn. The first five contributions, T\- 
T 5 , are products of the one-loop integrals defined in section |A.1| and are hence relatively 



straightforward to evaluate. We therefore simply define these contributions and present 
the corresponding results in terms of the integrals X, L, and Y. 



(N 2 - 1) 
16 



d 4 p sin 2 (|p i/ T) /sin 2 (ip M i?) 1 sm(^p^R) sin(±p M (i? - 2)) 



(2vr) 4 



D{p) V 



sin 2 (|p M ) 



(R + T)XL . 



+ 2 



sin 2 (^v) 



(A.20) 
(A.21) 



7 The five contributions in eq. ( |A.18| ) are denned explicitly in eqs. (3.6)-(3.10) of ref. 
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In eq. ( pOq ) 



d 4 k 



L 

2 



(2vr) 4 D(k) 

We repeat that we only keep terms which grow at least linearly with R and T. 



(A.22) 



T, 



(X 2 - 1) 



d 4 k 1 /sin 2 (|^T) ( lsin^i^i?) 
(2n)*D(k)\ sm 2 ak„) 



+ 77 



6 sin 2 (±£; M ) 



d 4 p Bin 2 (|p„T) sin 2 (|p M it!) 1 
sin 2 (|p M ) J 



(2vr) 4 £>(p) 



+ {( /U)J R)^(z/,T)} 



(AT 2 - 1) 



[XT-Y(T)] [Xi?-F(i?)] - 
(^"^r^o wom2 (^ 2 -l) 



96 



■[XR-Y(R)Y 



96 



[XT - F(T)] 2 . 



(A.23) 



(A.24) 



(X 2 - 1) 



12 

X 2 -l 
192 



d A p sin 2 (|p,T) sm 2 (l PlM R) 
(2tt) 4 D(p) sin 2 (ip M ) J 

{[XR-Y(R)} 2 + \XT-Y{T)} 2 } . 




+ {(n,R)~(v,T)} 



(A.25) 
(A.26) 



Ta 



2(X 2 -1) f d 4 k sm 2 (IKT) sin 2 (Ik^R) 



(2tt) 4 £>(jfc) 



sin 



2fl k 



X 



<i 4 p sin {\p^R) sin (t^T) 
(2vr) 4 D{p) sin 2 (i^) 



(X 2 - 1) 
24 



[XR-Y(R))(XT-Y(T)) . 



(A.27) 
(A.28) 



(X 2 - 1) / d 4 A; sin 2 (|fc MJ R) f d 4 p sm\±p v T) 



4 7 (27r) 4 J D(A;)sin 2 (|fcA J (2vr) 4 D(p) sin 2 (±p„ 
(X 2 - 1) 



16 



-(x J R-y(i?))(XT-r(T)) . 



(A.29) 
(A.30) 



The remaining three terms contain nested integrals and are considerably more compli- 
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cated to evaluate. The first of these is 



(N 2 



32 



d A k d A p 



(2vr) 4 (2vr) 4 D(p)D(k) 



sin(|(p + k)^) sin(|(p - k)^) - sin(|(p - k)^R) sin(±(p + k)^ 



x 



sin 2 (i(p + fc)„r) 
sin 2 (7^) sin 2 (|£;^ 



+ {(Ai,fl)<->(i/,T)} 



(A.31) 



In order to extract the behaviour of this integral at large R and T consider the following 
two-dimensional integral 



dky. dp^ 



1 



N 



2tt 2tt D(p)D(k) 4sin 2 (p M /2) 4sin 2 (fc /i /2) 4 sin 2 ((p + A;) M /2) 



(A.32) 



where 



iV 



4 sin 



(p + 



sin 



(P - &)/ 



4 sin 



(P - fc)^ 



sm 



(P + fc), 



It is convenient to use partial fractions to simplify I: 

I = h-{I 2 + h) + J 4 



(A.33) 
(A.34) 



where 


f dp^dk^ N 


Jl ~ 4A(p)A(k)J 
h = 1 ) 


2tt 2tt 4sin 2 (p M /2)4sin 2 (/^/2)4sin 2 ((p + k)j2) 
f dp^dk^ N 


4^(p)A(fc) J 


2tt 2tt 4 sin 2 (p M /2) (4sin 2 (^/2) + 2A(fc)) 4 sin 2 ((p + fc) M /2) 
^ dp^dk N 


/s ~ AA(p)A(k)J 


2tt 2tt (4sin 2 (^/2) + 2A(p)) 4sin 2 (/c M /2) 4sm 2 ((p + fc) /i /2) 
f dp^dkf, N 


h ~ 4A(p)A(k)J 


2vr 2tt 4sin 2 ((p+ fc) M /2) 



X- 



(4sin 2 (p M /2) + 2A{p)) (4sin 2 (fc M /2) + 2A(fc)) 



and A{k) are defined in eq. ( |A.2j) , where the sum is over the three components of p 
and k other than the /i-component. 

can be evaluated by changing integration variables to z = exp(ip M ) and uo = 
exp(i(p + k)fj,) and using Cauchy's contour integration theory. In terms of these complex 
variables 



AT = [(1 + z )(uj r - 1)(1 - uo/z) -{uo- l){z R + 1)(1 - {uo/z) R )] 2 . (A.35) 
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In this way we obtain 



h 

h 

h 
h 



R{R - 1) Y2 

4 

—XY(R) - ^X 2 + *XL 
4 v ; 4 4 

/ 2 



Note that 



N 



2tt 2tt 64sin J (p^/2) sin^fc^) sin 2 ((p + k)j2) 



R(R-V 



(A.36) 



Since we are only interested in identifying the terms which grow at least linearly with R 
or T, we can replace sin 2 (|(j> + k)„T) in eq. ( |A.31| ) with the factor |, thus obtaining 



N 2 



16 



3i? 



R 



X 2 XY(R) XL + RW 



+ R^T . 



(A.37) 



T, 



The second nested integral is 
N 2 -l ff d 4 k d 4 p 1 



6 

sm 2 (\p v T) 



(2tt) 4 (2tt) 4 D{p)D(k) 
sin(lp^R) sm^k^R) sin(|(p + k)^R) -y 
sin(|p M ) sin(|^) sin(|(p + k)^) 



} + {(/i, J R)^(z/,T)}J . (A.38 



Using the fact that 

d Pfl f dk„ sm( Pll R/2) sm(k^R/2) sin((p + k) ll R/2) 



2tt J 2tt sin(p M /2) sin(A; M /2) sin((p + k)j2) 



R 



(A.39) 



we readily obtain 



48 



X 2 (i? + T) 



(A.40) 
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The final contribution to Wi is 
(iV 2 -l) ff d 4 k d 4 p 



6 



(2tt) 4 (2tt) 4 D(p)D(k) 



sin 2 \\p v T)'- 



sin(|p M i?) rsin(|p At (i2 - 2)) cos((p + k)^) 



sind^) sin 2 (|£; p 



sin(|(p + k)^R - 2)) cos(|A; At (i? + 1) +p A1 ) sin(±^(i? - 2)) cos(|p M - fc 



sin(|(p + fc) M ) sin 2 (|/^) sin(ip A1 ) sin(±fc M ) sin(|(p + fc), 

sin(ifc M (i? - 2)) cos(±Qo + k)^R + - p M ) 



+ 



sin(|(p+ sin 2 (±/^) 
sin(i(p + fc) M (i2 - 2)) cos(iA; M (i2 + 2) + 



sin(±& M ) sin(|p M ) sin(±(p + 
sin(|A; M (i2 - 2)) cos(|(j9 - k)^ - ^ -p M )" 



+ {( yU ,i?)^(z/,T)} 



sin(ip M ) sin 2 (|/c M ) 

for which, by using partial fractions as in the evaluation of T$, we find 
„ (N 2 ~ 1) 



(A.41) 



48 



-X [-XR 2 -XR + 5Y(R)R + ALR] +R^T . 



(A.42) 



Summing up the terms T\ to T§ to obtain Wi we find 



(N 2 - 1) 



i=i 



48 



[XR-Y{R)){XT -Y{T)) 



(N 2 - 1) 
192 



{{XR - F(i?)) 2 + (XT - Y{T)f) 



(N 2 - 1) 

+ - - {-R 2 X 2 + 3RX 2 + 2XY(R)R - 2RXL + VIEW + R T} . (A.43) 

192 

W n : W/j is the "abelian part of the expansion of order g 4 " , and is defined in eq. (3.8) 
of ref. ||. It is very straightforward to evaluate: 



W 



n 



(2^-3)^-1) , 



6N 2 

(2N 2 -3)(N 2 - 1) 



96iV 2 

where W 2 has been defined in eq. (|1S|). 



[X(R + T) -Y(R) -Y{T)f 



(A.44) 
(A.45) 
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Wvp: Wyp is part of the contribution to the diagrams containing gluonic contributions 
to the vacuum polarisation, and is defined in eq.(3.9) of ref. ||. It can be written in the 
form: 



W 



(N*-l) (R + T) f d 3 p ( W^{p-p, = Q)-W^{p = Q) 

A 2 (p) 



VV N 8 J (2tt) 3 

where there is no implied sum over fi, and 



(A.46) 



N 



d A k 



(2n) 4 D(p + k)D(k) 



cos(^))(4 - C 4 (2p + k)) + sin 2 (^) (l + cos( Pl )) ] . 

(A.47) 



Although 11^ (p = 0) = 0, it is nevertheless convenient to subtract it (in the form of 
eq. ( |A.47| ) with p = 0) in the numerator of the integrand in eq. (|A.46|) and to write Wyp in 
terms of the integrals V\ and Vi defined in eqs. ( |A.7| ) and ( |A.9| ). This subtraction introduces 
non-zero contributions to each of V\ and V2, rendering these integrals convergent, but of 
course these contributions cancel in Wyp. In this way we obtain 



W VP = ^-^Wt) (V, + 2V 2 ) + ^^(R+T) - ^ 1 X . 



(iV 2 -l) /n , ^ /7L 5 



(A.48) 



Wyp: Wyp, the remaining contribution to the diagrams containing gluonic contributions 
to the vacuum polarisation M is defined in eq. (3.10) of ref. B and is simple to evaluate: 



Wrp = ^-3f- 1] HR + T)X - Y{R) - Y (T)] . 



96 A 2 

A. 2.1 Total Gluonic Contribution to W 4 : 



(A.49) 



Summing up the contributions from eqs. ( |A.43| ), ( |A.45[ ), ( |A.48| ) and (A.49), and keeping 
only those terms which grow at least linearly with R and T, we obtain for the total gluonic 
contribution to W±. 



W f uon = - ^_.}^ X\R + T) 2 + 



{N 2 



32N 2 



1QN 2 



-X(R + T)(Y(R) + Y(T)) 



fjV 2 - 1) 

+ - - {3X 2 - 2XL + 12W] (R + T) 



(A.50) 
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A. 3 Evaluation of the Fermionic Contribution to W 4 (R, T) 

In this section we evaluate the diagrams containing fermionic contributions to the vacuum 
polarisation. This contribution is also given by eq. flOED but with 11^ -> IT^ mion , where 
n^ mion is the fermionic contribution to the vacuum polarisation graph: 



11^ [p, p, - 0) - — J — s(kjS(q) ' (A ' 51) 



where q = p + k, S(k) is defined in eq. ( |A.11|) and Zw,sw(p, k) are defined in eqs. ( |A.13 ) 



and (|A.16|) respectively. Nf is the number of active light-quark flavours and we distinguish 
between the terms from the Wilson Fermion action (labelled by w) and the additional 
terms in the improved action with coefficient cgw (labelled by sw). Thus the fermionic 
contribution to W4 is 

fermion = ^ ~^ Nf (R + T) (V W + V SW ) . (A.52) 
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